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Abstract 

We study the approach to equihbrium for a scalar field which is coupled to a large 
thermal bath. Our analysis of the initial value problem is based on Kadanoff-Baym 
equations which are shown to be equivalent to a stochastic Langevin equation. The 
interaction with the thermal bath generates a temperature-dependent spectral den- 
sity, either through decay and inverse decay processes or via Landau damping. In 
equilibrium, energy density and pressure are determined by the Bose-Einstein dis- 
tribution function evaluated at a complex quasi-particle pole. The time evolution 
of the statistical propagator is compared with solutions of the Boltzmann equations 
for particles as well as quasi-particles. The dependence on initial conditions and the 
range of validity of the Boltzmann approximation are determined. 



1 Introduction 



The current standard model of cosmology explains many features of our universe as the re- 
sult of out-of-equilibrium processes during its very early high-temperature phase (cf. [1,2]). 
This includes the matter-antimatter asymmetry, i.e. the origin of matter, the production 
of dark matter, the formation of light elements and the decoupling of photons leading to 
the cosmic microwave background. 

Many nonequilibrium processes in the early universe can be treated in the canonical 
way by means of Boltzmann equations (cf. [1]) with sufficient accuracy. In some cases, 
however, quantum effects play a crucial role. This applies in particular to baryogenesis, 
the generation of the matter-antimatter asymmetry. Here the CP asymmetry, which leads 
to the baryon asymmetry, is the result of a quantum interference. It is therefore important 
to go beyond the classical Boltzmann equations and to treat the entire baryogenesis process 
quantum mechanically. 

An attractive baryogenesis scenario is leptogenesis [3,4], where a quantitative under- 
standing of the baryon asymmetry in terms of neutrino properties has been achieved [5]. 
In leptogenesis the out-of-equilibrium dynamics of a heavy Majorana neutrino, which is 
coupled to a large thermal bath of standard model particles, is the origin of the baryon 
asymmetry. Given the simplicity of this process, a full quantum mechanical treatment 
may be possible and some progress in this direction has already been made during the past 
years [6-8]. One important application is the study of flavor effects [9]. 

The treatment of nonequilibrium processes in quantum field theory is usually based 
either on Kadanoff-Baym equations and the Schwinger-Keldysh formalism [10-13] or on 
stochastic Langevin equations [14-17]. Both methods have been applied to various pro- 
cesses in particle physics and cosmology, including also electroweak baryogenesis [18]. In 
this paper we examine the connection between both approaches, which has been also con- 
sidered in [19]. As we shall see, the Kadanoff-Baym equations and the Langevin equation 
are, in fact, equivalent for the case of a large thermal bath where backreaction effects can 
be neglected. 

Boltzmann equations are first-order differential equations for number densities, which 
are local in time. They represent a valuable approximation for nonequilibrium processes 
in a dilute, weakly coupled gas. However, when the interactions between the quanta 
of the thermal plasma are strong, which is certainly the case in the presence of non- 
Abelian gauge interactions, the validity of the Boltzmann approximation is questionable. 
Correspondingly, the notion of number density becomes ambiguous, although several useful 
definitions have been suggested [13,17]. 

In this paper we study the approach to equilibrium for a scalar field which is coupled 
to a thermal bath with many degrees of freedom such that backreaction effects can be 
neglected. We shall focus on the description of this nonequilibrium process in terms of 
Green's functions rather than number densities. This is analogous to studies of preheating 
after infiation based on the statistical propagator [13,20]. As we shall see, the Kadanoff- 
Baym equations and the Langevin equation lead to identical results. 

Knowing the exact solution of the initial value problem for the Green's function of 



2 



the scalar field, we can systematically study the conditions for the validity of ordinary 
Boltzmann equations as well as Boltzmann equations for quasi-particles. At large times 
the scalar field reaches equilibrium. As we shall see, this state does not correspond to a gas 
of quasi-particles. There is an additional thermal 'vacuum' contribution which in principle 
can even lead to a negative pressure of low-momentum modes. The general solution of the 
Green's function also allows us to study the dependence of the equilibration on the initial 
conditions. This is an important problem in leptogenesis, because the baryon asymmetry 
can only be predicted in terms of neutrino properties when there is no dependence on the 
initial conditions [21]. 

To illustrate our results we consider a toy model of three scalars [17,22,23], one being 
much heavier than the other two. Two particles are in thermal equilibrium whereas the 
third one slowly approaches thermal equilibrium starting from zero initial abundance. Due 
to the interaction with the thermal bath this particle has a non-trivial spectral density, 
approximately described by a 'thermal mass' and a 'thermal width'. These are generated 
either by decays and inverse decays or by a process similar to Landau damping. Some 
aspects of this model have previously been studied based on the time evolution of a number 
density [17]. 

The paper is organized as follows. In Section 2 we define the various Green's functions 
in the Schwinger-Keldysh formalism and present a brief derivation of the Kadanoff-Baym 
equations. The theoretical framework leading to the Langevin equation is discussed in 
Section 3, following [17]. Section 4 deals with the solutions of the Kadanoff-Baym equa- 
tions. Thermal equilibrium and the quasi-particle picture are discussed in Section 5, and 
a sytematic comparison with Boltzmann equations is made in Section 6. The results are 
illustrated for a thermal bath of scalars in Section 7. A brief summary and outlook is given 
in Section 8. Various properties of the spectral function are discussed in the Appendix. 

2 The Schwinger-Keldysh formahsm 

Let us consider the nonequilibrium dynamics of a scalar field. In the Schwinger-Keldysh 
formalism the basic quantity is the Green's function defined on a contour C in the complex 
x°-plane (cf. Figure 1), 

Ac(xi,X2) =^cK,x°)A>(xi,X2)+^c(x^,x?)A<(xi,X2) . (2.1) 

The ^-functions enforce path ordering along the contour C, and A> and A< are the cor- 
relation functions 

A>(xi, X2) = ($(xi)$(x2)) = Tr(p$(xi)$(x2)) , (2.2) 
A<(xi, X2) = (<I>(x2)$(xi)) = Tr(p$(x2)$(xi)) , (2.3) 

where p is the density matrix of the system at some initial time U. 

We consider the case that the field $ is coupled to a thermal bath described by a 
self-energy 11. The Green's function Ac then satisfies the Schwinger-Dyson equation 

(□1 + m2)Ac(a;i, X2) + / d^x'nc(xi, x')Ac(x', X2) = -i5c{xi - X2) , (2.4) 

Jc 
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Figure 1: Path in the complex time plane for nonequilibrium Green's functions. 

where Di = {d"^ /dx\). Like the Green's function, also the self-energy can be decomposed 
as 

nc(xi,X2) =^cK,a;°)n>(xi,X2) + ^c(x^,a;?)n<(xi,X2) . (2.5) 

In the Schwinger-Dyson equation the time coordinates of Ac and lie can be on the 
upper or lower branch of the contour C, which we denote by the subscripts '+' and ' — ', 
respectively. Obviously, one has 

A_+(xi, X2) = A>(xi, X2) , A+_(xi, X2) = A<(xi, X2) , (2.6) 
n_+(xi,X2) = lP'{xi,X2) , n+_(xi,X2) = n<(xi,X2) , (2.7) 

whereas A_|_+, n++ and A , 11 are causal and anti-causal Green functions, respectively. 

From the Schwinger-Dyson equation (12. 4p one obtains for the correlation functions A< and 
A>, 

{Di + m^)A<{xi,X2) = y"rfV(-n++(a;i,a;')A<(a;',a;2) + n<(a;i,x')A„_(x',X2)) , (2.8) 

(□1 +m2)A>(xi,X2) = y"c/V(-n>(xi,x')A++(x',X2) + n__(xi,x')A>(x',X2)) , (2.9) 

where the relative sign in the integrands is due to the anti-causal time ordering on the 
lower branch of C. 

It is convenient to also introduce retarded and advanced Green functions, 

A«(Xi,X2) =^(tl-t2)(A>(Xi,X2)-A<(xi,X2)) (2.10) 
= ^(tl-t2)([0(Xl),0(x2)]) 
= A++(xi,X2) - A+„(xi,X2) 

= A_+(xi, X2) - A__(xi, X2) , 
A^(xi,X2) = -e{t2-ti){A>ixuX2)-A<{xuX2)) (2.11) 

= -^(t2-tl)([0(Xl),0(x2)]) 
= A++(xi,X2) - A_+(xi,X2) 
= A+_(xi,X2) - A„„(xi,X2) , 

n^(Xi,X2) = ^(tl -t2)(n>(Xi,X2) -n<(xi,x2)) 
= n++(xi, 0:2) - n+„(a;i, X2) 

= n_+(a;i, 0:2) - n__(a;i, X2) , (2.12) 
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n^(xi, = -d{t2 - ti)(n>(xi, X2) - n<(xi, X2)) 

= n++(a;i, X2) - n_+(a;i, X2) 

= n+_(a;i, X2) - n__(xi, X2) . (2.13) 

From Eqs. ( 12.8P and (12.91) one obtains the Kadanoff-Baym equations for the correlation 
functions A> and A<, 

{Di+m^)A>{xi,X2) = - J d^x' {n>{xi,x')A^{x',X2) + U^ixi,x')A>{x',X2)) , (2.14) 

-y"rfV (n<(a;i,a;')A^(x',a;2) + n^(xi,x')A<(x',a;2)) . (2.15) 



(□1 + m^)A<(xi,X2) = 

We now define the real symmetric and antisymmetric correlation functions 

A+(a;i, X2) = ^{{Hxi), $(0:2)}) , (2.16) 

A-(a;i,X2) =?([<f(xi),$(x2)]) , (2.17) 

and self-energies 

n+(xi,x2) = -^(n>(xi,x2) + n<(xi,x2)) , (2.18) 

n-(a;i,a;2) =n>(a;i,a;2)-n<(a;i,a;2) , (2.19) 

which also determine the retarded and advanced self-energies, 

U^{xi,X2) = 9{ti - t2)U-{xi,X2) , U^{xi,X2) = -9{t2-ti)U-{xi,X2) . (2.20) 

Adding and subtracting the Kadanoff-Baym equations (]2.141) and (12.15P , one obtains from 
Eqs. fl2.101) - (12.13p and (12.16p - (12.191) an homogeneous equation for A^ and an inhomoge- 
neous equation for A+, 

(□i +m2)A-(xi,X2) = - J d^^' dt'U-{xi,x')A-{x',X2) , (2.21) 
(□i +m2)A+(a:i,a:2) = - j d^^' j^' dt'U-{xi,x')A+{x',X2) 

+ j d^^' dt'U+{xux')A-{x',X2) . (2.22) 



We shall refer to these as equations as the first and second Kadanoff-Baym equation. A~ 
and A"*" are known as spectral function and statistical propagator (cf. [13]). Together they 
determine the path ordered Green's function, 



Ac{xi,X2) = A+(a;i,X2) - -signc(x'|' - xl)A~{xi,X2) . (2.23) 
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Figure 2: Path in the complex time plane for thermal Green's functions. 

A~ carries information about the spectrum of the system and A+ is related to occupation 
numbers of different modes. 

Using microcausality and the canonical quantization condition for a real scalar field, 



[^{x,)Mx2)]\t,=t, = [<i'(xi),$(x2)]|i,=t, = , (2.24) 
[$(xi), ^(xs)] \t,=t, = - X2) , (2.25) 



one obtains from the definitions (12.161) and (I2.17p 

A-(xi,X2)|ti=j2 =0 , (2.26) 

(9t,A"(xi,X2)|ti=t2 = -dt^A~{xi,X2)\t^=t2 = S{xi -X2) , (2.27) 

dtA2^'ixuX2)\t,=t,=0 . (2.28) 



In the following we shall restrict ourselves to systems with spatial translational invari- 
ance. In this case all two-point functions only depend on the difference of spatial coor- 
dinates, xi — X2, and it is convenient to perform a Fourier transformation. The Green's 
functions A^{ti,t2) satisfy the two Kadanoff-Baym equations 

rh 

{dl+ul)A-{h,t2) + / rft'n-(ti,t')A-(t',t2) = , (2.29) 
{dl+u;l)A+{h,t2)+ dt'U-{h,t')A+{t',t2) = dt'U^{tut')A-{t',t2) , (2.30) 
where u"^ = +m^. The initial conditions fl2.26l) - (l2.28p for the spectral function become 



A^{h,t2)\t,=t,=0 , (2.31) 
dt,A-{h,t2)\t,=t2 = -dt,A~{t,,t2)\t,=t2 = 1 , (2.32) 
dtA2^^iti,t2)\t,=t2 = . (2.33) 



For Green's functions in thermal equilibrium the density matrix in Eqs. (12.21) . (12. 3p 
is Peq = exp(— /9if), where H is the Hamiltonian of the system, and (3 = is the 
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inverse temperature. Time coordinates of Green functions now lie on the contour shown 
in Figure 2, and one has invariance under time translations so that two-point functions 
only depend on the time difference ti — ^2- After a Fourier transformation, one obtains the 
KMS relations [12] for Green's functions and self-energies, 

- coth ( — 
2 V 2 



A+M = --coth(^^jA-M, (2.34) 
n+(c.) = -%othf^')n-(^). (2.35) 



The Kadanoff-Baym equations describe the dynamics of a arbitrary nonequilibrium 
system. Depending on the self-energy and the initial conditions, the solutions will generally 
be complicated. An enormous simplification is achieved for a large medium such that the 
backreaction of the field $ can be neglected. Furthermore, we assume that the medium is 
in thermal equilibrium and, therefore, the self energy of $ is time-translation invariant, 

nq(tl,t2) = nq(ti-t2) . (2.36) 

In this case also the spectral function is time-translation invariant, as shown in Ap- 
pendix A.l. With these simplifications, the Kadanoff-Baym equations become 



{dl + -t2) = - r dm-^{t, - t')A-{t' - t,) , (2.37) 

Ju 

- I ' dt'T{-{U - t')A+(t', h) . (2.38) 
Jti 

These equations will be solved in Section 4 for general initial conditions. 



3 Stochastic Langevin equation 

Nonequilibrium processes can also be studied by means of Langevin equations which 
describe the evolution of the field itself rather than the evolution of Green's functions 
(cf. [14-17]). Below we sketch a brief derivation of the Langevin equation describing a 
scalar field $ coupled to a large thermal bath with bosonic and fermionic fields x, follow- 
ing the discussion in [17]. We assume that the coupling is of the form and neglect 
the backreaction of $ on the thermal bath, which makes the problem solvable. 
The starting point is the nonequilibrium generating functional [13, 17] 
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where the subscript 'in' stands for the initial condition. The action of the fields $ and x 
is given by 



5[$±, X, J±] = / d^x + g^+0[x+] + J+$+ 

Ju 



where is the Lagrangian of a free massive field 

C 



1 1 

m <P 



= 2 

and Pin stands for the matrix elements of the initial density matrix, 

Pin($+;$'rj = ($|p|$'). 



(3.2) 



(3.3) 



(3.4) 



The field $ lives on the Keldysh contour C shown in Figure [U $±(x) is the field with the 
time argument on the "forward" (C+) and "backward" (C_) part of this contour, respec- 
tively, satisfying the boundary conditions 



$+(t„x) = $+(x) , <l>_(t„x) = <l>r^(x) . 



(3.5) 



The fields x ^ire assumed to be in thermal equilibrium, corresponding to the contour Cfs 
(Figure EJ, which is possible since the backreaction of $ on the thermal bath is neglected. 
In the following we shall choose as initial time tj = 0. 

It is convenient to perform a change of variables in the functional integral (]3.ip , 



R{x) = $+(x) - <l>_(x) . 



(3.6) 
(3.7) 



We are interested in the two-point function of \E', which couples to the source term J 
— Integrating out the fields R and x one finds [17], 
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Z[J] = jD^u^Dnu,W{^in,ni^) J V^V^V[^]e' 

Jo 

here the measure V[C,] is given by 

(-1 /"OO /"OO 

2 Jo Jo ^^'^'l(^)nq(^-^0-'e-q(t' 



(3.^ 



(3.9) 



and 'Cq(^) is a stochastic noise. The Fourier transform ^q(t) in (j3.8l) satisfies the initial 
conditions 

^q(O) = ^q,in , ^q,in(0) = n^,in • (3.10) 
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The function >V(\E'in; Ttm) is a functional Wigner transform of the initial density matrix, 

M^in; n,^) = I Divine- /'^^--W^-Wp,, (^^,, + - . (3.11) 

For a pure vacuum state p is a product of the two delta functions (5(\E'in) and 5(7rin). 

In order to obtain two-point correlators of the field \1/ one has to solve the classical 
stochastic Langevin equation, 

{d^ + ool) ^^{t) + r dt'U^it - t')^^{t') = Ut) , (3.12) 
Jo 

with the initial conditions ()3.10l) . Since the backreaction of the field $ is neglected, the 
only relevant correlation functions are 

iUt)) = , (3.13) 
iUm^'it')) = -n+(t-t')5(q + q') . (3.14) 

The solution of the Langevin equation is conveniently expressed in terms of an auxiliary 
function /q(t) which is defined as solution of the homogeneous equation 

{d^ + ^l) m + r dt'U^it - t')Ut') = , (3.15) 

^0 

with the initial conditions 

/q(0) = , /q(0) = 1 . (3.16) 

One easily verifies that the solution of the Langevin equation is then given by 

*q(i) = *q,m/q(t) + ^^Mt) + dt'f^{t - t')^^{t') . (3.17) 

Correlation functions of the scalar field can now be obtained by calculating the expec- 
tation values 

(vl/q^(ti)...vl/q„(t„)) , (3.18) 

which involve the correlation functions of the stochastic noise and also an average over the 
initial conditions. For the simplest case, the two-point function, one has 

{'^ci{h)^^'{t2))=g^{ti,t2)6{q + q')=g^{t2,ti)6{q + q!) . (3.19) 

From the Langevin equation (I3.12P one easily derives an analogous equation for the two- 
point function, 

{d^ + a;^) (v^q(ti)vl/q,(i2)) + r dt'U-it, - t')(^q(t')^q'(^2)) (3.20) 

^0 

= (eq(tl)^q'(t2)) (3.21) 



9 



= 5(q + q') / ' dt'U^ih - t')Ut' - h) , (3.22) 
Jo 

which implies 

{dl + uj''^g^{hM)+ r dm-^{h-t')g^{t'M) (3.23) 

Jo 



= ['dt'U+{t,-t')f^{t'-t,) . (3.24) 

A solution of this equation can be directly obtained from the solution of the Langevin 
equation f l3.12p . In the case where the initial field and its time derivative vanish, 

= = , (3.25) 

the relevant averages for the two-point function are 

(^q,m^q,m) =5(q+q')«q , (3-26) 
(^q,in^q',in) =5(q+q')/5q , (3-27) 
(^q,in^q,in) =5(q+q')7q • (3.28) 

Using the solution (13.171) and the correlations f l3.14p one obtains the two-point function 

^q(tl, h) = «q/q(tl)/q(t2) + lj{h)f{t2) (3.29) 
+ /3q (/q(tl)/q(t2) + /q (^1 ) /q(t2 ) ) (3.30) 

+ r dt' r dt"Uh - t')U+{t' - t")Ut" - h) . (3.31) 

Jo Jo 

In the following section we shall see that the auxiliary function /q(t) and the two-point cor- 
relation function 5'q(ti, ^2) are precisely the spectral function and the statistical propagator 
of the field $, respectively. 

4 Solving the Kadanoff-Baym equations 
4.1 The equation for the spectral function 

As proven in Appendix A.l, the spectral function is time translation invariant, i.e., it only 
depends on the time difference y = ti — ^2- Hence, the first Kadanoff-Baym equation (12.37p 
takes the form 

{di + ^l)Kiy)+ rdy'U-{y-y')A-{y') = 0. (4.1) 
^0 

This equation can be solved by performing a Laplace transformation, 
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A-(.) = / dye~^yA-{y) , (4.2) 
Jo 
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with 



for which one obtains after a straightforward calculation 

roo poo 

nj(s)=/ e-'yU^{y)dy= e-^^n-(|/)dy = fl-(.) . (4.4) 
Jo Jo 

According to (14.31) . the general solution of (14. ip depends on two parameters, the values of 
A~ and dyA~ at y = 0. Using the inverse Laplace transform one finds 

A,fe) ^ aA,(0) + A-(0)a,) £ I^^J^^I^ . (4.5) 

Here Cb is the Bromwich contour (see Figure 3): The part parallel to the imaginary axis 
is chosen such that all singularities of the integrand are to its left; the second part is the 
semicircle at infinity which closes the contour at Re(s) < 0. Since the integrand of (14. 5p has 
singularities only on the imaginary axis, the second part can be deformed to run parallel 
to the imaginary axis as well: Cr — » _|_ r 

The spectral function A~(|/) satisfies the boundary conditions (I2.3ip and (I2.32p . which 
implies 

A-{y) = / — . (4.6) 

This result can be further simplified by making use of the analytic properties of the self- 
energy n~(s). On the real axis n~(s) is real, while on the parts of the contour which are 
parallel to the imaginary axis one has 

n-(icu ± e) = RenJ(cu) ± ilmnj(cj) , (4.7) 

with 

ImnJ(..) = i [U^iu + le) - Yiliuj - le)) . (4.8) 
Hence, the expression (14. 6p takes the form 

J-oo -^TT 

where the spectral function Pq(c<j) is given in terms of real and imaginary part of the 
self-energy Yi^iuj)^ 

-2Imn^(cj) + 2uje 

Note that Imnq(ci;) and Re 0^(0;) are odd and even functions, respectively, which implies 
that Aq(|/) is real. Further properties of this solution are discussed in Appendix A. Let us 
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Figure 3: Bromwich contour 

recall that the expression ( 14.101) is obtained after neglecting the backreaction of the field 
$ on the thermal bath. This is the reason why the self-energy and the spectral function 
are time translation invariant. 

The self-energy U^i^uj), and consequently the spectral function Pq(co'), are divergent 
and have to be renormalized. This can be done by the usual mass and wave function 
renormalization at zero temperature. In fl4.10p is replaced by i^q(o) = ''^o + Q^i where 
mo is the bare mass of the field $. The difference between bare and renormalized mass 
squared is determined by requiring that at zero temperature the spectral function has a 
pole at cjq = + q^, 

- <o) - RenJ(^q)|r=o = . (4.11) 
Expanding the self-energy around around Uq, a further divergence can be absorbed in a 
wave function renormalization constant, 

ReU^iu) = ReU^{u^)\T=o + (l - Z~') [u^ - cu^) + ReflJ(u;) , (4.12) 
where Renq(ti;) is the finite part and 



2co'q du 



(4.13) 



The spectral function (14.101) now takes the form 



Pq(^) = Z- ^ . (4.14) 



(u^ -LU^- ZReU^iuj)^ + {ZlmU^iu) + cue) 



2 



Introducing the renormalized field operator = y/Z^, one obtains the renormalized 
spectral function p'^{uj) = Zp^{uj) in terms of the renormalized self-energy Il^''^{uj) = 
^nj(a;), 

-2Imn^'''(a;) + 2u;e 
(u^ -uol- Ren^'^cu)) + (imn^'^cu) + cue' 
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The divergencies of spectral function and statistical propagator can be removed in the same 
way by mass and wave function renormalization at zero temperature. In the following we 
shall drop the superscript 'r' to keep the notation simple. 

The spectral function describes a quasi-particle resonance at finite temperature with 

nl-iul-ReU^i%) = 0, nl\T=o = iol, (4.16) 

and decay width 

^ -i-lmnj(l]q) . (4.17) 

For simplicity, we have neglected the effect of Imllq on the quasi-particle energy. The 
correction 6flq = O^T"^) is evaluated in Section 6. 

In a free theory Im 11^(0;) = 0, and (I4.15P is a representation of the 5-function. The 
spectral function (14.91) then oscillates without damping, i.e., there are no dissipative effects. 
Dissipation arises either from $ decays and inverse decays or, similar to Landau damp- 
ing, from scattering processes with particles in the plasma. Which of these mechanisms 
dominates the dissipative effects and therefore the equilibration process depends on the 
position of the quasi-particle pole relative to the masses of particles in the thermal bath. 
A specific example will be discussed in Section 7. For small width the spectral function 
is well approximated by the Breit-Wigner function. The relevant formulae are collected in 
Appendix A. 5. 



4.2 Solution for the statistical propagator 

We are now ready to solve the second Kadanoff-Baym equation (12.381) for the statistical 
propagator, which for initial time = is given by 

idl+ujl)A+ih,t2) + / rft'n-(ti -tOA+(t',t2) = C(ti,t2) , (4.18) 

Jo 

with ^ 

C(ii, t2)= r rft'n+(ti - t') Aq (f - h) . (4.19) 

One easily verifies that the solution can be expressed as 

A+(ti,t2)=A+(ti,t2)+ r dt'A-{h-t')C{t'M) , (4.20) 
where A+(ti,t2) satisfies the homogeneous equation 

{dl + ujl)kX{hM) + / ' dt'Ii-{U - t')A+it',h) = . (4.21) 

The homogeneous equation is identical to (14.11) . with t2 playing the role of a parameter. 
We can therefore read off the general solution from (14.51) . 

A+(ti,t2)=Aq(t2)A-(tO+5q(t2)Aq(tl) . (4.22) 
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Using the symmetry A^(ti,t2) = Aq(t2,ii), one obtains 

Aq(t2)Aq(tl) + i?q(t2)A-(ti) = A^{t,)A-{t2) + 5q(tl)A-(t2) • (4.23) 

Together with the boundary conditions (i23TD - (l233l ). Aq(0) = Aq(0) = and A-(0) = 1, 
this implies 

A^{t) = A^{0)A^{t) + i?q(0)A^ (t) , B^{t) = iq(0)A-(t) + i?q(0)A^ (t) . (4.24) 

Inserting A^it) and i?q(t) in fl4.23p and using the symmetry of A^ (ti, ^2), one finds -Bq(O) = 
Aq(0). The initial state of the system is therefore characterized by three constants, which 
can be chosen as 

Kin = A+(ti,t2)|t,=i,=0 = Aq(0) , (4.25) 
A+.„ = dt,A+ih,t2)\t,=t,=0 = 9i,A+(ti,t2)|i.=t,=0 = i?q(0) = iq(0) , (4.26) 

AJ.„ = dtA.Kitut2)\t,=t.=o = B^{0) . (4.27) 

From Eqs. K22h . and the initial conditions ([L5l)-(0) we now obtain the 

full solution for the statistical propagator, 

A+(ti,t2) = A+i„Aq(ti)Aq(t2)+A+i„Aq(tl)Aq(t2) 

+ A+,„ (A-(ti)A, (t2) + Aq (ti)A^(t2)) 

+ KmeJtut2) , (4.28) 

where 

A^„(^i,^2)= f'dt' r dt"A-{t,-t')U^{t' -nA~{t" -t,) . (4.29) 
Jo Jo 

This contribution to the statistical propagator, which is independent of the initial condi- 
tions, is often referred to as memory integral. It can be expressed in the form 

A:„(^i, ^2) = - r (*-*^)H;(ti, u;)n^{h, u;)U+{u) , (4.30) 

where [17] 

n^{t,u)= I rfre-^-^Aq(r) . (4.31) 
Jo 

The expression f l4.30p will be the basis of our numerical analysis in Section 7. 

5 Thermal equilibrium and quasi-particles 

Let us now verify that the solution (14.281) for the statistical propagator approaches thermal 
equilibrium at late times. This means that the quantity 
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which becomes a Fourier transform for t ^ oo, satisfies the KMS condition asymptotically, 

A+(oo,^) = -^cothf^') . (5.2) 



For late times only the memory integral is relevant, since A^(t) and A^{t) fall off 
exponentially for t ^ l/F. One then obtains 

A+(oo,a;) = A+„(oo,a;) = -\n^icK^,uj)\'U+{u;) . (5.3) 

The quantity 7Yq(oo,ct;) is the Laplace transform of the spectral function, 



POO 

n^{oo,u)= / cire-'('"-*^)"Aq(r) 
Jo 



A-{tuj + e) 



s=iui+e 
1 



-ujI- Ren^(o;) - ilmU^{u) 



(5.4) 



which yields 



YHJoo, u 



|2 



''^ (cu2 -ul- Ren«(c^))2 + (Imn^(a;))2 

Pq(^) 



5.5) 



2Imn^(cj) ■ 

Inserting this expression into (15. 3p . using the KMS condition for the self-energy and (IA.30I) . 

nq(a;) =2zImnJ(o;) , 

one obtains (cf. M.l liinl )). 



At{oo,uj) = — coth 



\ o / (,.,2 



Imn«(cj) 



2 ; (cu2 - - Ren^(cj))2 + (Imn^(cu))2 
%othf^')A^(^). (5.6) 



2 V 2 ^ 

Hence, our solution for the statistical propagator indeed fulfills the KMS condition (12.34p 
in the limit t ^ oo, which proves that the system reaches thermal equilibrium. For a 
specific example the approach to equilibrium will be studied numerically in Section 7. 

It is instructive to evaluate the statistical propagator in thermal equilibrium at equal 
times, i.e., y = h — t2 = 0, 
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For a free field one has 

Pq(u;) = 27rsign(u;)(5(u;^ - u^) , (5.8) 
which yields the well know result 

A+|,=o = ^Q + riBK)) , (5.9) 
with the temperature dependent Bose-Einstein distribution function 

^bK) = . (5.10) 

Generically, the interaction with the thermal bath changes the energy cUq of a free 
particle to a temperature dependent complex energy fiq which appears as a pole of the 
spectral function Pq{uj) and the integrand of (15.71) . The spectral function then has two 
poles in the upper plane, fiq and — ^^q, which are determined by the condition 

n^-[iul + u^{n^)Y' = 0. (5.11) 

Assuming that the integral can be closed in the upper half-plane, one obtains for the 
statistical propagator in equilibrium^ 

A+|,=o = Re(^^Q + nB(^^q)) I . (5.12) 

Compared to (15. Sp . the Bose-Einstein distribution function has been replaced by the com- 
plex distribution function nB(^q)- 

At high temperatures, where <^ 1, the Bose-Einstein distribution has a well-known 
infrared divergence, 

ubM ~ > 1 . (5.13) 

pU^q 

For quasi-particles, where c<jq is replaced by fiq = fiq + this divergence is cut off by 

the finite width, 

which remains finite even if the real part fiq vanishes. 

Comparison of equations (15. 9p and (I5.12p suggests that in thermal equilibrium the $ 
particles may form a gas of quasi-particles. This question can be clarified by evaluating 
energy density and pressure of the $ particles. Since the expectation value of $ vanishes, 
one obtains from the energy momentum tensoiEl 

T^. = 9,$ - 7/^,L (5.15) 



^Here we restrict ourselves to the case where there are no additional poles. 
We use the convention diag = (1, — 1, — 1, — 1). 
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for the contribution of a mode with momentum q to energy density and pressure, 

= (Too)lq = \{^' + iy^f + m'^')U , (5.16) 

= = (^(V$)^ + - (V<|.)^ - m^<|.2))|, . (5.17) 

This yields for the energy density 

ec,(00) = \ {dtA.+^D Aj(tl,t2)|,^=,^_ 

= ^K + ^^)^Q + ^B(f^.)) , (5.18) 



and for the pressure 



3 

where, for simplicity, we have neglected the quasi-particle width. 

In summary, the energy momentum tensor in thermal equilibrium can be expressed as 
sum of a quasi-particle gas contribution and a temperature dependent 'vacuum' term, 

(T^.) U = n,u, (ef + ) - %uPf + r/^.<^^ . (5.20) 
Here u^' = (1,0) is the 4-velocity of the thermal bath, and 

«f = (5 + «B(Si<,)) , (6.21) 

Energy density and pressure of the quasi-particle gas agree with the corresponding expres- 
sions for a free gas, with the energy tUq of a free particle replaced by the quasi-particle 
energy f2q. The 'vacuum contribution' k^^^ vanishes for f2q = cUq. For large thermal 
effects, i.e. r2q ^ LJq or f2q <ti c^^q, the equation of state differs significantly from the one 
of a free gas. Note that for < cj^, the pressure can even become negative! 



6 Comparison with Boltzmann equations 

The time evolution of nonequilibrium systems is usually studied by means of Boltzmann 
equations for particle number densities. However, this notion does not have a well defined 
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physical meaning in a nonequilibrium process. For a dilute, weakly coupled gas the number 
density of 'free particles' may be a good approximation, and in some cases the the effect of a 
medium can be taken into account by considering quasi-particles. In general, however, one 
has to study the time evolution of Green's functions, in particular if quantum interferences 
are important. 

In order to determine the range of validity of the Boltzmann approximation, we shall 
consider in this section the time evolution of an observable, the energy density. The exact 
expression can be obtained from the statistical propagator, and approximations are given 
by solutions of Boltzmann equations. In this way, the description of the nonequilibrium 
process by means of Green's functions on the one hand, and Boltzmann equations on the 
other hand, can be directly compared, based on the same observable. 

Consider the Boltzmann equation for a dilute gas of $ particles. The competition 
between a gain and a loss term determines the change of the particle number density [22], 

dtn^it) = (1 + nq(t))7< - nq(t)7> , (6.1) 

where production and decay rates satisfy the KMS relation and are obtained from the 
self-energy of the field 

7q =e-^"''7q =riBK)7q, (6.2) 
7q = • (6.3) 

Using these relations, the Boltzmann equation (16.11) can be written in the form 

dtUqit) = -7q(nq(t) - njiiuj^)) , (6.4) 

with the obvious solution 

ridit) = riBiLU^) + (nq(0) - nB(c^q)) e'^^' . (6.5) 

For comparison with the Kadanoff-Baym equations we now consider instead of the 
number density the energy density of a mode with momentum q, normalized to the energy 
of a single quantum, 

eq(t) ^ !aW = 1 ^^(i) . (6.6) 
The deviation from the equilibrium density. 



^q(i) = + 5^q(t) , (6.7) 

with 

eq(oo) = 6j- = i + nB(a;q) , (6.8) 
satisfies the differential equation 

(dt + 7q)5eq(t) = . (6.9) 
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The modification of the spectral function in a thermal bath (cf. (14.10p ) suggests to 
replace the equilibrium value and the evolution equation for the energy density by the 
expressions 

eq(i) = + Se^it) , = 1 + n^{n^) , (6.10) 

and 

{dt + T^)6e^{t) = , (6.11) 
where the quasi-particle width is given by (cf. Appendix A. 5) 

Imn^(fi„) , 1 9 p , 

= V = i-^a;R<M|„,. (6.12) 

As long as the interaction of the field $ with the thermal bath can be treated perturbatively, 
the difference between solutions of the two Boltzmann equations for particles and quasi- 
particles, respectively, should be small. When the quasi-particle width becomes large, 
however, the use of first-order differential equations, which are local in time, becomes 
clearly questionable. 

As discussed in Section 5, the exact time dependence of the energy density can be 
directly obtained from the statistical propagator, 

'-"^^^ = + ^q) K(tiMt,=t,=t ' 

which satisfies the Kadanoff-Baym equation (12.150 . 

{dl+ul)A+{h,t2)+ [\t'U-(ti-t')A+{t',t2)= dt'U+{h-t')A-{t'-t2) . (6.13) 

For large times, t ^ l/Fq, the dependence on the initial values at = can be neglected, 
and one obtains 

/oo 
dt'{U^ih-t')A^it',h)+tU^ih-t')A^it'~h))=0. (6.14) 
-oo 

Changing time variables. 



^^h+t2^ y = ti-t2, A+(t;2/) = A+(ti,t2) , (6.15) 



and expanding. 



t' - = A+ it; t' - + ^5tA+ (t; t' - t2) + . . . , (6.16) 
one finds for the Fourier transforms with respect to the time differences, 

■^l-^^^t-uJ^ + ul^ A+(t;^) 

19 
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Using the relations (jATzl) - (jXaTl) . one obtains from the real and the imaginary part 
of this complex equation two equations for the real quantity Aq(t,u;), 



(^^dl - uj^ + A+(t,^) = -RenJ(a;)A+(t,.;) + n+MlmA^(t, 



iaimn^(cj)9A+(t,cj) 



a;^A+(t,^) = ImnJ(a;)A+(t,a;) + n+HReA^(t,a;) 

iSRen^H 9A+(t,cj) 
^2 di ' 

where the dots indicate neglected higher-order terms. 

Consider now an expansion around the equilibrium solution, 



(6.19) 



A+(t,a7) = A+H + 5A+(t,^) . (6.20) 

From equation fl6.19p one reads off 

ImnJ(cu)A+(cu) + n+(cu)ReA^(cu) = , (6.21) 

which is satisfied because of ([Al30l) . (1X221) and the KMS conditions ([2:341) and (l235|l . 
The first equation (I6.18P yields for the equilibrium solution, 

{u' -ul- Renj(c.)) A+(.;) = -n+(^) Im A^ (a;) . (6.22) 

In the zero- width limit, this equation is fulfilled for 



= fiq = ^ujI + Re n^(f^q) . (6.23) 
The finite width leads to a correction, 

uj = Vl^ + 5Vl^. (6.24) 
Expanding (|6.22l) in (5f2q, one obtains to leading order 

2nq5fiqA+(fiq) + n+(l]q) Im A^(fiq) = , (6.25) 

which implies 

„ r,(n,)imA;;(n,) 

We can use the free spectral function, 

A„ (cj) = 27rsign(cu)(5(tu2 - fi^) ^ ^6 27) 
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to evaluate Im Aq(f2q) to leading order in Fq, 



toA>,)^4p/-^..'^^. (6.28) 



Using fICTll . flXMl) . fl4l(l and (TOTll we finally obtain 

5fiq = • (6-29) 

Hence, for Fq ^ fiq, the leading term in the derivative expansion indeed implies oj = fiq. 
If finite width effects are not negligible, however, off-shell effects become important and 
the derivative expansion becomes unreliable. 

Inserting u = Qq^ in the first-order differential equation ()6.19l) . one obtains for the 
departure from equilibrium of the statistical propagator, 

- 2?i;£^«-'l«.) I - *Aj(i;n,) = o . (6.30) 

Hence, (9^(5 A+(t; fiq) = 0{lmll^), and to this order Eq. ( 16.181) is also satisfied. 
We can now evaluate the energy density 

^~qW = ^ {dtA. + ^D K(^^Mt.=t,=t 

= du(^^d^ + u;' + u;l^A^{t;u;), (6.31) 

which approaches the equilibrium value 

^^'°°'^^"' = ^(2+"B(nj) . (6.32) 

As already discussed in the previous section, the true equilibrium value of the energy 
density does not correspond to a gas of quasi-particles, 

^ ■ (6-33) 

From Eqs. (I6.30p . (lA.4ip and (IA.42I) one obtains for the deviation from the equilibrium 
value, 

{dt + Fq) eq(t) = , (6.34) 

which is identical to the Boltzmann equation (16. lip for quasi-particles. 

In summary, we have obtained the following conditions under which the approach to 
equilibrium can be described by Boltzmann equations. For a dilute, weakly coupled gas 
the ordinary Boltzmann equation for a number density is sufficient, which approaches the 
Bose-Einstein distribution for a gas of free particles. When interactions with a thermal 
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Figure 4: Spectral function Pq(c<j) for q = 0; case (a) with masses mi = m2 = 0.2m and 
temperatures Ti = 0.1m, T2 = 0.2m, T3 = 0.5m. 

bath significantly change the spectral function, a Boltzmann equation for quasi-particles 
describes the approach to equilibrium as long as the quasi-particle width can be neglected. 
However, the equilibrium value of the energy density is different from the one for a gas of 
quasi-particles. Finally, when the width cannot be neglected and off-shell effects become 
significant, a linear evolution equation of 'Boltzmann type', which is local in time, is 
no longer adequate. Instead, the dynamics is non-Iocal in time, and one has to solve 
Kadanoff-Baym equations. 



7 A thermal bath of scalars 

So far we have performed a very general analysis, and the only approximation has been to 
neglect the backreaction of the field $ on the thermal bath. Furthermore, we have restricted 
our discussion to the case that $ is linearly coupled to the bath via an interaction term 
g^O{x) (cf- (13. 2p ). In general, x represents an arbitrary number of bosonic or fermionic 
fields with arbitrary couplings including gauge interactions. In order to illustrate the results 
of the previous sections, we now consider a toy model (cf. [17,22]), where the quanta of 
two massive scalar fields represent the thermal bath. The full Lagrangian is given by 

C = ^d.^d'^^ - ^m2$2 + ^ ^^^,x^^''Xi - Irnjxl) + 9^X1X2 + C^int • (7.1) 

Note that the coupling g has the dimension of mass. In the following we shall neglect 
self-interaction of the x fields and use free thermal propagators for simplicity. 

We consider two cases: (a) m ^ mi,m2 and (b) m2 ^ m,mi. In the first case, 
dissipation is dominated by $ decays and inverse decays, $ ^ X1X2, whereas in the second 
one X2 decays and inverse decays, X2 ^ ^Xi, are most important. 
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Figure 5: Spectral function Pq^uj) for q = 0; case (b) with masses mi = m, m2 
temperatures Ti = m, T2 = 2m, T3 = 5m. 
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Figure 7: Real part of the self-energy Il^{uj) for q = 0; case (a) with masses nii = m2 = 
0.2m and temperatures Ti = 0.5m (solid) and T2 = m (dashed). 



In both cases the imaginary part of the self-energy is known analytically [17,23]. The 
relevant formulae are collected in Appendix IA.3[ For m ^ mi, 7712, the decay width of $ 
at zero temperature is given by 



g^2 



^ m. (7.2) 



IGtt \m 

To illustrate thermal effects we shall use a rather large coupling which corresponds to 
T/m = 0.1. 

The spectral function Pq(co') (cf. ()4.15l) ) is shown in Figures [4] and [5] for the two mass 
patterns (a) and (b), respectively. In case (a), 11^ has an imaginary part at zero tem- 
perature. The width is large, and already at small temperatures the quasi-particle profile 
becomes broad. On the contrary, in case (b) the zero-temperature width is zero and the 
finite-temperature width is small. Hence, the quasi-particle profile becomes broad only at 
much larger temperatures. The spectral function A~{y) is the Fourier transform of 2pq(cij). 
As Figure [6] illustrates, it approximately represents a damped oscillation with frequency 
flq and damping rate Fq. 

It is interesting that thermal corrections can increase or decrease the particle mass 
m. Whether the quasi-particle peak moves to the right or to the left depends on the 
position of the zero-temperature pole relative to the branch cuts, and it also depends on 
the temperature. This can be seen by considering the real part of the self-energy, which is 
displayed for two different temperatures in Fig. [7] for case (a). For the smaller temperature 
one has Renq^Q(m) > 0, whereas Renq^Q(m) < holds for the larger temperature, which 
corresponds to a shift of the particle mass to the right and to the left, respectively. 
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Figure 8: Statistical propagator Aq(ti,t2) for q = 0; case (b) with masses mi = m, 
1712 = 5m and T = 10m. 

The statistical propagator A^{ti,t2) depends on the initial conditions. The most gen- 
eral gaussian initial density matrix has five free parameters (cf. [13]). We consider the 
simplest case of a free field density matrix and vanishing mean values $ and $, which 
implies for each momentum mode, 

<f q,in = , (7.3) 
<^'q,in = , (7.4) 

Kin = A+(ti,t2)|t,=fe=o = ^{l + ^<i) ' (7-5) 
Kin = 5*iA+(ti,t2)|ti=t.=o = dt,A+ituh)\t,=t,=o = , (7.6) 
Kin = dtA2Kitut2)\t,=t,=o = uJ^(^l + n^^ . (7.7) 

The initial state of the system is now characterized by only one parameter which 
corresponds to an initial number density for a free field. 

The general solution ( 14.28^ for the statistical propagator (ti, ^2) is shown in Figure [3 
For fixed t = (ti+t2)/2 one sees damped oscillations iny = ti — ^2- The amplitude increases 
with increasing time t, as illustrated by Figure M For fixed y = ti — t2 one observes the 
approach to equilibrium with increasing t = {ti + t2)/2. For large times the departure from 
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Figure 9: Statistical propagator Aq(ti,t2) as function of y = ti — t2 for q = 0; case (b) 
with mi = m, 1712 = 5m, T = 10m and three values of t = (ti + t2)/2: mt = 15 (dashed 
line), mt = 20 (dotted-dashed), mt = 60 (solid). 



A^(ii,^2) 




Figure 10: Statistical propagator A'^{ti,t2) as function oft = {ti + t2)/2 for y = 0, q = 0; 
case (b) with masses mi = m, m2 = 5m and T = 10m. 
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Figure 11: Statistical propagator A'^(ti,t2) with q = as function of t = (ti + 12)/2 for 
y = and different initial conditions; case (b) with masses q = 0, mi = m, m2 = 5m and 
T = 10m. 

equilibrium is described by a first-order differential equation, and it decreases exponentially. 
At small times the evolution is governed by a second-order differential equation, which leads 
to the oscillations visible in Figure [TOl The independence of the equilibrium solution from 
the initial conditions is illustrated by Figure [TTl The memory of the initial conditions is 
lost at times t > 1/T. 

Finally, it is important to recall that the equilibrium value of the energy differs from 
the one obtained in the Boltzmann approximation. This is illustrated in Figure [l2] where 
the different contributions to the energy are compared as functions of temperature. For 
the chosen parameters the particle and quasi-particle energies are indistinguishable. The 
'vacuum contribution' is positive, which means that the total energy is larger than the 
particle/quasi-particle one. The reason is that for the chosen parameters thermal correc- 
tions decrease the particle mass. For other parameter choices the 'vacuum contribution' 
can have opposite sign. 

8 Conclusions and outlook 

We have studied the approach to equilibrium for a real scalar field coupled to a large 
thermal bath. We have computed the exact two-point functions, the spectral function and 
the statistical propagator, for arbitrary initial conditions. This is possible for a thermal 
bath with many degrees of freedom such that the backreaction of the scalar field can be 
neglected. 
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Figure 12: Energy density eq = Cq/cjq as function of temperature for q = 0; case (b) with 
masses mi = m, m2 = 5m: total energy density (solid), particle and quasi-particle energy 
densities (dotted), and 'vacuum' energy density (dashed). 

The self-energy representing the thermal bath is time-translation invariant. We have 
shown that this is also the case for the spectral function, whereas the statistical propagator 
depends on two time coordinates, ti and t2, and also the time ti where the initial conditions 
are specified. 

We have obtained the two-point functions by solving the Kadanoff-Baym equations, 
which turned out to be equivalent to solving a stochastic Langevin equation. As expected, 
the relaxation time is determined by the imaginary part of the self-energy, i.e., a 'quasi- 
particle width' r. For t > l/F, the statistical propagator becomes independent of the 
initial conditions. It is then given by a memory integral which depends on the real and 
imaginary part of the self-energy. 

As long as thermal corrections are small, the approach to equilibrium is well described 
by the ordinary Boltzmann equation, which is a local, first-order diff'erential equation in 
time for the particle number density. However, in the case of large thermal corrections 
the notion of number density becomes ambiguous, and it is important to consider the 
second-order Kadanoff-Baym equations for the two-point functions rather than a Boltz- 
mann equation. Still, as long as the quasi-particle decay width is small compared to the 
quasi-particle energy, a Boltzmann equation for quasi-particles describes the approach to 
equilibrium to good approximation. For large decay width the dynamics becomes nonlocal 
in time and the Boltzmann approximation breaks down. 

It is interesting to study the contribution of the thermalized scalar field to energy 
density and pressure. For a free field these observables are determined by the Bose-Einstein 
distribution function. Interaction with the thermal bath can significantly modify energy 
density and pressure, and therefore the equation of state. The Bose-Einstein distribution 
as function of the complex quasi-particle pole is now the relevant quantity. Energy density 
and pressure differ from the expressions for a free gas of quasi-particles by a temperature- 
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dependent 'vacuum term' which can become important at high temperatures. 

We have illustrated these results for a toy model where the thermal bath consists of 
two massive scalar fields. We have considered two cases where equilibration takes place 
either via decays and inverse decays or via 'Landau damping'. In general, one has to 
study the Kadanoff-Baym equation for the statistical propagator which depends on two 
time coordinates as well as initial conditions. However, for large times, t ^ l/F, the time 
evolution is well described by the Boltzmann equation for quasi-particles. 

Our analysis has been motivated by the need of a full quantum mechanical description 
of leptogenesis. To achieve this, one has to consider correlation functions rather than 
number densities, although for parts of the calculation the use of Boltzmann equations 
will be sufficient. The heavy Majorana neutrino is very weakly coupled to the thermal 
bath. Hence, thermal corrections to its mass and width are small, and its approach to 
equilibrium is well described by Boltzmann equations. However, to study the dependence 
of the final baryon asymmetry on initial conditions it may be necessary to consider 
the statistical propagator, since leptogenesis takes place at ts ~ l/F. Furthermore, for 
lepton and Higgs fields, which have strong gauge interactions, finite-width eff'ects can 
be important. At present it is unclear how accurately the leptogenesis process can be 
described based on a quasi-particle picture for the standard model particles which form 
the thermal bath. These questions are currently under investigation [24]. 
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A The spectral function 
A.l Time-translation invariance 

In this section we shall prove that the most general solution of the first Kadanoff'-Baym 
equation is time-translation invariant. The starting point is Eq. (12 .29^ with the boundary 
conditions (12.3ip - (|2.33l) . Performing the change of variables ti = t + y/2, t2 = t — y/2, 
Eq. (ICTl ) becomes 

Qs,^ + dtdy + dl+ ujI^ A^it; y) + dy'U^iy - y')A^it'; y') = , (A.l) 

where t' = t — {y — y')/2 and = + rn? . Note that and 11" only depend on |q| 
because of rotational invariance. Both functions are antisymmetric in y. The boundary 
conditions f l23TD - (l233l ) read 

A-(t;0)=0, (A.2) 
aiA-(t;0)=0, (A.3) 
a,A-(t;y)|,=o = l , (A.4) 

^^t'-^j) Aq(t;y)|,=o = 0. (A.5) 

The condition (IA.2I) is automatically fulfilled because of the antisymmetry in y. 

To prove that A~ is time-translation invariant we now perform an expansion in powers 

of n-, 

oo 

^q=E^f^ AW = 0(nW). (A.6) 

n=0 

For n = one has 

(\^dl + d,dy + dl + ujI^ a W (t; = . (A.7) 
Using the antisymmetry of A^ in y, one obtains 

dtdyA^^\t-y)=Q, (A.8) 

which has the general solution 

A?)=<)(t) + &fHy)- (A.9) 

Every solution of Eqs. (1A.7I1 and (lA.SP satisfies the boundary condition (lA.Sp . The condi- 
tion ( 1A.3I) implies 

dtA^^\t-Q)=dta^^\t) = Q . (A.IO) 
Hence, Oq "* is constant and Aq'* only depends on y. Eq. (1A.7I ) now becomes 

(^.' + ^3 AjHy) = , (A.ll) 
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which has the antisymmetric solution 

A^^\y) = c^^hmico^y) . (A.12) 
For n one can use the recurrence relation 

(^\d^ + dtdy + dl + A("+^)(t; y) + dy'U-^{y - y')A(-\y') = . (A. 13) 

Using the antisymmetry of 11" and A~ in y, one again finds 

dtdyA^^+'\t;y) = 0. (A.14) 

Repeating the same steps as for Aq yields the result that also Aq"'*'^'' is independent of t. 
We conclude that the spectral function is the antisymmetric solution of the equation 

(9j + c.2)A^(y)+ rdy'Il~^iy-y')A^iy') = 0, (A.15) 
Jo 

with the boundary condition 

dyA^{y)\y=o = I . (A.16) 

A. 2 Conventions for propagators and self-energies 

In thermal equilibrium the retarded and advanced propagators and self-energies only 
depend on the time difference y = ti — t-2. In the following we list several relations between 
their Fourier transforms, which are used in the different sections. In principle, these 
relations are all well know, but their specific form depends on the chosen conventions. All 
relations are not affected by the three-dimensional Fourier transform. We therefore drop 
the argument q or x. 

Propagators: 





= -A-(^) 




(A.17) 


A+(cu)* 


= A+H, 




(A.18) 


^\^) 


= -A"(cu)- 


^ du;'A-{u') 


(A.19) 


A\u) 


^ A / 

= — A~(u 
2 ^ 


^ ^rdu'A-{j) 
y„oo uj' -u 


(A.20) 


ReA^{u) 


= -ReA^( 


^) = ^A-H , 


(A.21) 


lmA^{uj) 


= lmA^{u] 


7-00 2vri a;' - a; ' 


(A.22) 


A^i-u;) 


= A^H. 




(A.23) 
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(A.24) 



Self-energies: 



n-(w)* = -n-(cj) , (A.25) 

U+{u)* = U+{u) , (A.26) 

n» = -Vm + ^/?^^^, (a.27) 



2 J 2TXi LO' — LO 

2 J Zm uj' — uj 

Reli\u) = ReU\u;) = V f^^^^ , (A.29) 

J 2m u — u 

Imn^(u;) = -lmn^(cu) = ^n-(cu) , (A.30) 

n^(-u;) = n^H . (A.31) 
A. 3 Scalar field model 

The interaction with the thermal bath changes the spectral function of a free scalar particle, 

Pq(cu) = 27rsign(cu)5(cu^ - cUq) , (A.32) 
to the expression (14.151) which depends on real and imaginary part of the self-energy, 



-2ImnJ(cu) + 2u;e 
^''^'^^ " _ ^2 _ Ren^(cu)]2 + [Imn^(w) + ue]^ " ^"^'^^^ 



We have computed the imaginary part of the self-energy in the scalar field model defined 
in Section d assuming free thermal propagators for the fields xi and X2- The result agrees 
with [17]. One obtains {q = (cjq, q)): 

- Imnj(c.) = ao{q) + af{q) + a^^{q) . (A.34) 
Here ctq is the zero-temperature contribution due to the decay process $ X1X2, 



o"o(g) =77 — ^sign(u;)e(g^ - (mi + ma^^^ 



X 



((g2)2_2g2(^2^^2^ + (^2_^2)2^2 ^ ^^35^ 



a^'^'* is the finite-temperature contribution from this process, 

= ^g^|q|^ sign(a;)e(g^ - (mi + ma)^) 
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1 — e~^'^+ 

X ]+im,^m,)] , (A.36) 



and o"?^(q) is the finite-temperature contribution from processes Xi ~^ 



b g"^ 
^^^^^^^ ^16^^^^^^"*^^'*®'^*^"^^ ~ "^^^^ ~ 

( (\- f.-ti\^-\\ 

X In + ("^1 ^ ^2) ) , (A.37) 



1 _ e-/^l^+l 



where we have used the abbreviations 



The real part of the self-energy can be computed using the dispersion relation relation, 

1 ,Imn^(cu') 
Ren^(cu) = -V / dJ . (A.39) 




Figure 13: Poles and cuts of the spectral function p{uj) for q = at T = 0: (a) m > mi+m-^, 
and (b) m < mi+ 1712- 

Based on these expressions we can discuss the analytic structure of the spectral function. 
For a free field Pq{uj) is given by (1A.32I) which has two poles at u; = ±u!g in the complex u- 
plane. The interaction of $ with Xi and X2 does not modify these poles for m < mi +m2, 
where $ is stable at zero temperature. In addition there are branch cuts at the two- 
particle thresholds \uj\ > cuthi = v^q^ + {rrii + 777,2)^ (see Fig. [T3b). They correspond to 
virtual decays and inverse decays, $ ^ XX- In the case m > mi + m2 these processes can 
happen on-shell since m > c^thi, and $ becomes unstable. Now the spectral function has 
four poles in the complex ci;-plane, whose real parts lie in the region of the branch cuts (see 
Fig. [T3h). The imaginary parts of the poles correspond to the decay width of $. 

^Note that we disagree with the discussion in [22] which imphes the additional factor 6(|mf — mj] — g^). 
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(a) 



+ Ren^ 



(6) 



Figure 14: Poles and cuts of the spectral function p{uj) for q = at T ^ 0: (a) m > mi+m2, 
and (b) m < mi + m2. 



The analytic structure of the spectral function at finite temperature is displayed in 
Fig. [m The position of cuthi is shifted due to thermal corrections from Rellq. Fur- 
thermore, a new branch cut appears in the region where ai, ^ 0, i.e. for \uj\ < a;th2 = 
+ (mi — 7712)^. This is due to processes x 4'X ^^nd corresponds to Landau damping 
of quasi-particles in the plasma. If the real part of the poles falls into the regions of one of 
the branch cuts, i.e. < Ci;th2 or > ciJthi, they acquire an imaginary part which corre- 
sponds to the quasi-particle decay width (see Fig.[T4|). Qualitatively, this analytic structure 
is typical for interacting quantum field theories at finite temperature. In general, the spec- 
tral function can have additional singular contributions for m < |a;| < cuthi corresponding 
to bound states. At finite temperature they are also dressed to quasi-particles. 



A. 4 Breit-Wigner approximation 

In the regime of couplings and temperatures where | Imnq(f2q)| <C f2q, so that the quasi- 
particle picture holds, one can approximate the spectral function Pq{uj) by a Breit-Wigner 
function. From the expression (14. 15^ one easily obtains 



where Fq is the quasi-particle width 
with 



sign(t^)rc 



2fiq(| 



UJ\ 



l-p2 
4^ q 



Imng(l]q) 



1 (9Ren^(cj) 



2a 



duj 



(A.40) 



(A.41) 



(A.42) 
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Contrary to the exact spectral function (I4.15p . the Breit-Wigner approximation (1A.40I1 
has no branch cuts. The integrals over uj are dominated by the regions around the quasi- 
particle poles where the two functions are very similar. For the Fourier transform, the 
spectral function in real time, one obtains 



A-(,)^Z,^^^e-^^*/^ (A.43) 
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